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Abstract 

Additional variables (also often called "hidden variables") are sometimes added to standard 
quantum mechanics in order to remove its indeterminism or "incompletness," and to make the 
measurement process look more classical. Here we discuss a case in which an additional variable 
arises almost spontaneously from the quantum formalism: the emergence of relative phase between 
two highly populated Fock state Bose-Einstein condensates. The model simulated here involves the 
interference of two Bose condensates, one with all up spins, and the other with down spins, along a 
z-axis. With the clouds overlapping, we consider the results of measuring spins in a transverse plane 
(the general direction is studied in an appendix). The determination of the previously "hidden" 
phase becomes progressively more definite as additional measurements are made. We also provide 
an analysis of a recent and closely related experiment. 
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I. INTRODUCTION 



Several laboratories have made use of Bose condensates with spins or pseudo-spins to 
perform very interesting experiments. For example, in JILA experiments a mixture of 
two hyperfme states playing the role of pseudo-spin- 1/2 particles was found to segregate by 
species, thereby exhibiting spin waves. Another "spin" experiment j^j] involves overlapping 
two Bose condensates, one with spin up and the other with spin down and observing the 
appearance of a spontaneous transverse spin polarization. In this paper we consider a 
simulation that involves the interference of two clouds of Bose gases, one with up spins 
and the other down spins along some z-axis. With the clouds overlapping, we consider 
what happens in measuring spins along a transverse direction at a set of azimuthal angles, 
detecting whether a spin is found up or down along each measurement angle used. (The case 
of a general measurement direction is studied in Appendix A.) In looking at a simulation 
of this process we note an interesting element involved in the process: the relation of the 
successive measurements to the appearance of a so-called "hidden" variable. 

The introduction of additional variables to the standard formalism of quantum mechan- 
ics is not a new idea; it dates back almost to the appearance of this theory [3J, with, for 

n n 

instance, the early work of L. de Broglie 5] and later D. Bohm on "hidden variables" p. The 
main motivation was to restore determinism by reproducing the statistical results of quan- 
tum mechanics with classical averages over additional variables. In 1935, Einstein, Podolsky 
and Rosen (EPR) showed that there exists an even stronger argument to complete quan- 
tum mechanics, without referring to determinism: the standard form of quantum mechanics 
does not satisfy local realism, and should be completed by additional "elements of reality" 
to restore it U Thirty years later J. Bell, with a famous theorem 0, jj, extended the EPR 
argument and showed that, not only the formalism, but also the predictions of quantum 
mechanics are sometimes incompatible with local realism. In other words, adding variables 
to quantum mechanics is not sufficient to restore local realism in all cases: there exist some 
situations where the evolution of additional variables has to be explicitly non-local. This 
remarkable result stimulated several generations of careful experiments in order to decide 
whether or not quantum mechanics still gives correct predictions, even in these surprising 
"non-local" situations. An impressive body of evidence has now been accumulated in favor 
of the predictions of quantum mechanics, even if none of these experiments is ideal 10]; the 
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general consensus among physicists is that these surprising "non-local" predictions of quan- 
tum mechanics are indeed obeyed by Nature. This does not mean that additional variables 
should now be excluded from quantum mechanics! Actually, advocates of these variables 
argue that it is precisely one of their big merits to make the non-local character of quantum 
mechanics explicit |l3j|- For a general discussion of the consequences of the Bell theorem, 
see for instance Refs. 111! I2I 

In most cases, the introduction of additional variables into quantum mechanics leaves 
a large range of flexibility: one has to introduce whatever new variables seem appropriate 
(positions, momenta, fields, etc.) with adequate equation of evolution, chosen so that a 
statistical average over initial condition restores the usual predictions of quantum mechanics. 
There is nevertheless a case where an additional variable emerges almost spontaneously, in 
an unique way, from the quantum formalism: the spontaneous appearance of relative phase 
between two highly populated Fock states (Bose-Einstein condensates). The general physical 
phenomenon was discussed long ago by Anderson |l4(, with more emphasis on spontaneous 
symmetry breaking in phase transitions than additional variables in quantum mechanics. 
Later several authors used various formalisms, often borrowed from quantum optics, to 
give more detailed calculations I^J. The approach used in Ref. is slightly different 
since it shows how a simple conservation rule, the conservation of particle number, can 
naturally be expressed through an integral over the "conjugate" variable, the relative phase 
$ of the two states. The probability of any sequence of results then appears as a sum over 
this new variable, exactly as in theories with additional variables. Moreover, each time a 
measurement is performed, the state vector projection postulate provides a new initial state, 
which changes the $ distribution; it turns out that the change of this probability distribution 
can be obtained very easily for any sequence of measurements. In other words, one has access 
to the evolution of the distribution function of the additional variable under the effect of one 
or more successive quantum measurements. This is what we study in the present article; for 
the sake of simplicity, we limit ourselves to the limit of large occupation numbers, a case in 
which the Bell theorem does not predict any incompatibility between quantum mechanics 
and local realism, so that locality will not be an issue here. 

Additional variables in quantum mechanics are often called "hidden variables" for histor- 
ical reasons. Neverheless, J. Bell pointed out how inappropriate this name is ("Absurdly, 
these theories are known as hidden variable theories...." )|22|], pointing out that the name 
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would be more appropriate for the standard wave function of quantum mechanics. Indeed, 
in an interference experiment, for instance, these additional variables are not hidden but 
actually directly observed in the result of the individual experiments. On the other hand, 
the wave function or state vector can be reconstructed only indirectly by statistical analysis 
after many measurements. This is why we will tend to avoid the words "hidden variables" 
and rather speak of "additional variables" here. 



II. PHASE IN SPIN STATES 

Suppose we have particles with two internal states, either real spin 1/2 or pseudo-spins 
as in the case of two hyperfme states. If we have two stationary clouds of Bose particles 
with N + spin-up, and AL spin-down, particles along a z-axis, the initial state is 

\^) = \N + ,N_), (1) 

a Fock state in spin space. We want to measure the occurrences of a sequence of spin 
measurements in the transverse xy-pl&ne, at a series of azimuthal angles, <px, <p2, • • • <Pm, 
resulting in the sequence of results {rji} , either up (+) or down (— ) along the angles: 

7 7l = ±l; 77 2 = ±l;---7 ?m = ±l. (2) 

Here we present a simplified calculation of the probabilities, which com plet ely ignores orbital 
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We have a set of 



variables; a more precise calculation is given in Appendix A and Ref. 
angular momentum variables given in terms of the destruction operators a for particles up, 
and b for down, along a z-axis. The number operator is then 



n = a^a 



+ b\ (3) 



with the angular momentum variables 

a z = a) a — b^b, 
a x = a)b + tfa, 



a 



i{tfa-a)b). (4) 



These definitions should make clear the meaning of the axes of our xyz coordinate system 
and the term "transverse." The z-axis is defined by the two states (possibly hyperfme states) 



into which our particles are condensed in Eq.(l). Thus a z is diagonal in Fock states, while 
a x and a y are not. For example, while |1, 0) as an eigenstate of a z , an eigenstate of a x is the 
mixed state |1, 0) + |0, 1). 

The expectation value of the operator 



2N 
1 

2iV 



71 



f] (cos 4>ia x + sin 4>i(Jy)} 



n + rjU^a + e'^a^b 



(5) 



gives the probability of finding a spin with component rj along the angle 0, in the transverse 
plane. (To see the connection between this operator and the probability, see Appendix A or 
Ref. [21I ) Suppose one measures a single spin along a transverse axis at angle 0i starting in 
the state \§>) quoted above. The probability of finding spin with result 771 is easily seen to 
be 
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P 1 (</> 1 ) = {N + ,N-.\ — 



1 



\N + ,N.) = - 



(6) 



as expected. However, if immediately after, one measures a second particle along a different 
transverse axis 02 then a straightforward calculation gives the probability of the sequence 
{771,7/2} to be 



P 2 (0i,0 2 ) = (N + ,N_\p m (fo)p m (4>i)\N + ,N_) 



x 



1 +7?i?72y cos(0i 



(7) 



where x = 2yfN + N_/N and N = N + + AL. Having found the first particle up or down 
along (f)\ affects the result of the second measurement along 02. The two measurements are 
correlated by boson statistics. It is possible to write this last result in an instructive way as 



1 r 2 * ci$ 

2 ) = 7 / tt( 1 + cos (^i - $ )) ( X + XT l2 cos(0 2 - $)) 

4 JO ZTT 



(8) 



as one can verify by doing the integration. It is remarkable that this form of the probability 
for a sequence of m such measurements along transverse axes 0i, 02, • • -0m persists (2 1|: 

P - ~ / o~ II I 1 + ^ cos (^ - $ )] • ( 9 ) 

This result just quoted assumes all measurements are done in the transverse plane. How- 
ever, measurements could also be done of spin up and down along an arbitrary axis at angles 
(pi). The generalization of Ref. I21I to this case is described in Appendix A. The derivation 



there, when specified to all Oi = n/2 and uniform orbital variables, reduces to Eq. as 
expected. 

The form of Eq. (jUJ) is quite interesting. Suppose we had started out with a system 
polarized partially (at a fraction x) along the transverse direction $; then the resulting 
probability for finding a set of spins up or down along the set of angles {0j} would be Eq. 
but without the integration over $. Since we started out knowing the numbers of particles 
in longitudinal states, all transverse phases must be present equally in the initial state as 
represented by the integral over the transverse phase in Eq. Q. 

We can look at the expression of Eq. is a somewhat different way: The probability 
that the r\ in the mth spin measurement is ±1, after a sequence of results {r]i}, is 

P m {±)~ d$2 m ($)(l±xcos(0 m -$)) (10) 
Jo 

where the mth measurement is made only in the transverse plane at the angle m . In this 
we have 

m—1 

9m($) =11(1 + xrji cos(0 l - $)) . (11) 

i=l 

The function g m ($), as we will see by explicit simulation, peaks up, after a sufficiently 
large number m — 1 of measurements, at some value of phase, call it $ - Thus subsequent 
measurements will appear as if the the system of particles had, after just m — 1 measurements 
(with m <C N), been prepared with a set polarization phase angle. (Of course, if we repeated 
the experiment, starting again from the same initial state, a different phase would emerge 
and indeed it is random in a series of such experiments.) Because particle number and 
phase are conjugate variables, the original state of known particle number has, after m 
measurements, morphed into a state of relatively well-known phase, but with the number 
of particles up or down along z much less certain (although the total number or particles 
remains known). 

We start with = \N + , AL) . Our first measurement along some azimuthal angle 
produces the new state 

= ^{n + T)(e i(t, tfa + e- i ' t> a}b)) \N + ,NJ). (12) 

It is clear that this process produces a mixture of states \N + , iV_) , |iV + + 1, iV_ — 1) , and 
\N + — 1, AL + 1) , Each subsequent measurement produces a further mixing of number states 
so that number becomes less certain. Its conjugate, phase, becomes less uncertain. 
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As an alternative view one might think of the phase as having been there all along, but 
temporarily hidden from view, with experiments continually clarifying its value. In such a 
view one integrates probabilities over all possible values of an additional variable introduced 
to "complete" quantum mechanics, as we find in Eqs. (JOJ) or (jiup . The additional quantum 
variable is this "emerging" phase angle. 

In the rest of this section we confine our discussions to the simplified case of N + = AL = 
N/2, or x = 1. Let us now Fourier transform g m (Q) to write 



g m ($) = < + £ [< cos (<? $ ) + K sin(g$)] . (13) 

3=1 

In Eq. (jl()|) . only the q = 1 term contributes, so that 

h m \ 1 

(14) 



Pm{±) = ~ 



1 , f a T , K . 

1 ± cos <p m H sm . 



2aff 2a r 



With this result we can define two parameters that characterize the the results of the 
mth measurement. These parameters are <3> m and a m (or A m ) given by 



cos <£> = ; ; sin<P m = , =: sma m = A r , ■ 



a? 



■m i j ° 1J - L ^ m I j in — in „ r ,, 

(15) 

(Because P m is a probability, A m < 1 and we can write A m as a sine.) We obtain 

P m (±) = - [1 ± sin a m (cos 4> m cos $ m + sin m sin $ m )] . (16) 

Consider a single spin polarized as some space angle (/3, 7) and assume that a measurement 
if performed along a transverse direction <fi m . The probabity of finding a + result is given 
by the well-known expression 

|± (4>m I A 7) I 2 = 2 t 1 ±sin/5(cos0 m cos7 + sin m sin 7)] , (17) 

which has exactly the same form as Eq. (fTB^) if /3 = a m , 7 = $ m . Actually the two expressions 
are equal whenever the "spin" lies on a cone of directions around the measurement direction 
<f) m containing this particular direction. We show this "cone of equal probability" (CEP) in 
Fig. ^ with a spin at arbitrary (3, 7 on the cone. 

What this calculation shows is that two different points of view are possible. At each 
measurement step, all the effects of the previous measurements are contained in the distri- 
bution function g m (<&); in this first point of view, the system is described by a statistical 



mixtures of spins polarized in transverse directions. But one obtains the same probabilities 
for the next measurement by replacing this statistical mixture by a single pure spin state, 
fully polarized along directions (3 = a m , 7 = $ m (or any direction making the same angle 
with the direction of measurement). A second point of view is therefore possible, where the 
spin is fully polarized, but in a direction that is, in general, no longer in the transverse 
plane. For instance, if the distribution g m (&) has no $ dependence, the $ m dependence of 
the probabilities disappears, and this spin is polarized in a direction that is perpendicular to 
the direction of measurement (a m = 0). On the other hand, the $ m dependence is maximum 
when a m = n/2 so that the pure spin state lies in the transverse plane. 

III. THE AMHERST EXPERIMENT 

Ref. reports an experiment that is well- described by our analysis. In Ref. U the ex- 
periment is treated in terms of a two-component spinor, with each component representing 
one of the two condensates in a hyperfme state with a well-defined phase. Here we consider 
the experiment in terms of a Fock state having an up-spin-down-spin ratio the same as the 
experiment. The Amherst experiment combines two condensates each originating from a 
different point and claims that the result shows a "spontaneous transverse polarization." 
However, what they are able to detect by laser absorption is the existence of each of the lon- 
gitudinal components. So after mixing the two components they perform a 7r/2 tip in order 
to measure the size of the polarization and its phase. What they see then is an anticorrelated 
interference pattern in the two components. 

We assume a Fock state for the initial mixed-condensate function: \^) = \N + , AL) , with 
an arbitrary pre-established x, y, z coordinate system. We introduce a set of transverse 
coordinates it, w offset by an angle cfi relative to the original xy axes: 

it = cos (fix + sin <fiy (18) 
w = — sin (fix + cos (fiy (19) 

We want to measure the w-component of the spin in multiple simultaneous measurements 
made along the longitudinal (z) direction. To enable this we do a spin tip of tt/2 around 
the w axis. Then one measures spins along z. So the wave function analyzed is 

\W) = U \V) = U \N+, NJ) (20) 
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where 

U = e inaw/4 (21) 

and 

<7 W = — sin 4>a x + cos (j>a y (22) 

After applying U we make measurements along z, that is, we look for the probability 
given by 

m i 

p m = m & n 2N {n + Vi(Tz) u m (23) 

This is precisely equivalent to the relation 

m i 

Pm = (*l II 2V (» + *U*<r,U) |tf > (24) 
i=i 



But we have 



U^a z U = siii(f)ay + cos<pa x (25) 



so that 



P m = (N+, N-lH — in + rji (sin0a y + cos^)) \N + , AL> (26) 
i=i 



1 f 271 " 
2 



/ 2 "rf$ p m ($)(l + r/ m xcos(0-$)) (27) 

JO 



where the last line follows from the discussion of the previous section. We have seen that if 
we make simultaneous measurements over the whole cloud of particles g m ultimately peaks 
sharply at a particular phase $o , so that 

P m (v m ) ~ ^(l + r] m xcos(0-$ o )) (28) 

which is just the result given in Eq. (6) of Ref. U. What this shows is that the up and down 
measurements will be anticorrelated in phase. Ref. uses this anticorrelation in phase in the 
spatial interference patterns of the two detected clouds. That reference assumed two point 
sources of the two bosonic clouds released from traps. The discussion of Appendix A does 
indeed include a spatial phase; if we use 9{ = n/2 there, we find we find that a spatial phase 
£(r) should be added to the argument of the cosine in Eq. (J2EJ)) where £(r) = arg [iz (r)/w&(r)] 
and u a (r) and Ub(r) are the single particle wave functions describing the two condensates. An 
order of magnitude of the behavior of the interference fringes can be derived by assuming we 
are allowed to substitute simple one-dimensional spreading Gaussian wave packets for u a (r) 



and Ub{r). (Appendix A does not explicitly mention time dependence in the wave functions 
but that inclusion is straightforward.) When we do this we find that the time-dependent 
phase is of order £(r,t) = md 2 /ht, where d is the separation of the centers of the two wave 
packets, m the particle mass, and t is the time after the initially narrow packets have been 
spreading. This feature of time-dependent fringes is assumed in the Amherst experment as 



it was in Ref. 
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23l A complete derivation of this fringe result is given by Wallis et al. 
Our interpretation of the Amherst experiment shows the ambiguity of interpretation that 
can arise in such an experiment. The interpretation given by Ref. [l] is that the spontaneous 
polarization existed in the transverse plane before the ninety degree tip took place. But did 
it? Before the longitudinal measurements, did $o exist? One might consider it an additional 
variable that was already there, but needed to be brought out by experiment. Or one can 
say it did not exist; the wave function was a linear combination of all such phases and 
was collapsed by longitudinal measurement to the final value $o found in the end. Either 
interpretation works. 



IV. NUMERICAL SIMULATION 



The physics becomes clearer if we do a numerical simulation of the measurements. The 
first spin result is chosen randomly up or down along fa . To choose from the probability P m 
at each subsequent stage (m = 2, 3, • ■ ■) we simply pick a random number p from to 1. If 
p < P m {+), then 7] m = +1; if p > P m {+), then r} m = -1. 

The first set of experiments is performed with all measurements done at the same angle 
fa = 0. However, we find that to gain any more information about the emerging phase angle 
we must switch measurement angles. 



A. Experiment with constant direction of measurement. 

Our results are shown in the first figures for a run of 300 trials with all fa = fa = 0. There 
is no plot of the azimuthal angle $ m because generally in the case of a constant measurement 
angle, $ m will be equal to fa or fa + ir for every iteration. This is easy to understand 
physically by a symmetry argument: with a constant direction of measurement, there is 
no way to distinguish between spins polarized in two directions that are symmetricial with 
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respect to the direction of measurement in the transverse plane. With <\>\ = 0, (or consider 
the symmetry in terms of $' = $ — <pi), g m {§) must be a symmetric function around $ = 0; 
by Eqs. (|T3j) and f(15j) we see b™ = 0, sin$ m = 0, and so $ m = or n. For an alternative 
view of the same result, look at Eq. (|16J1 in the case of all <p m identically equal to <\>\. This 
equation is of the form 

P m (±) = - [1 ± sin a m cos(0! - $ m )] (29) 

But every cosine in Eq. ()l()j) is a function of $ — <pi, and so the integral is independent of 
0i, as one finds by changing integration variable to $ — <\>\. From Eq. (J2U|) this implies that 
$m — <pi — constant, which we show in Appendix B must be or it. 

Let m± be the number of spins found with 77 = ±1 along <\>\. We must always have just 
Pm{+) — |(1 =t sina m ). Since we anticipate physically that P m (+) — m + /{m + + m_) for 
large m, the plus sign occurs when m + > m_ and the minus sign in the opposite case. In 
this constant 4>i case, the variable A m = sin a m goes to some unpredictable random positive 
value as seen in Fig. EJ (As it does so, it has slightly curved trajectories followed by small 
jumps: the trajectories occur during a series of identical r] values; a jump happens when the 
opposite 1] eventually occurs. We explain this peculiar behavior in Appendix C.) Fig.Elshows 
#300 ($) versus $ (as well as the earlier and wider g 10 and gi 50 ) . It peaks up sharply at two 
locations, near ±$0 with $0 — 0.73, symmetrically situated around azimuthal measurement 
angle <p\ — 0; this is exactly what we should expect, as mentioned above, since, in keeping 
the same measurement angle, the probability is identical for spins symmetrically located to 
the sides of measurement angle. 

There is another piece of data from the run, the number of 77's of each sign. We find that 
there were m + = 262 up spins and m_ = 38 down spins in the experiment. When m + > m_ 
as we have here, we have approximately 

^(l + sina m ) = ^± (30) 
2 m 

with m = m + + m^. Solving we get sin a = 0.75, which is very close to what we get from the 
limiting value in Fig. El In Appendix C we show that one can do the g integral analytically 
for constant measurement angle, and the actual results for sin a is remarkably simple, namely 

I TYli j TTIj I 

sina m = (exact result for constant 6). (31) 

m + 1 
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Numerically this is 0.74. For large m the values of Eqs. (j3"Uj) and (|31|) are identical since the 
1 in Eq. (|3*T|) is then negligible. 

The above results of a m — > constant, and $ m = 0, seem to imply that the "hidden 
polarization direction," whose orientation we are trying to determine, is at zero azimuthal 
angle and either above or below the transverse plane at polar angle a m . This angle in turn is 
related to the angles ±$ of the peaks of g m (®) through the cone of equal probability (CEP) 
of Fig. The cone opening angle is $o and the polar angle of the cone is the complement 
a = 7r/2 — $ - A general position on the CEP around the measurement angle 0i = obeys 

P(±) = | ± (0|/?, 7 )| 2 = ^[l±sin/?cos 7 ] = ^[l±sin«] 

= i[l±cos$ ] (32) 



, 7T 

1 ± sin - - $ 



.2 

Thus all we can say about the hidden-spin angle so far is that sin j3 cos 7 = sin a = cos$ - 
We could have 7) = (a,0) , or (/3,j) = (7r/2,$ ) or somewhere else on the cone, but we 
cannot tell yet which it is. The distribution <?($) peaks at the two places where the CEP 
crosses the transverse plane. 

We can analytically show the double peaking property of <?($) in the case where all 
measurement angles 0, are equal, say, to zero. Ref. |l8| shows how this function has two 
equal sharp maxima in the range we consider. For completeness we repeat this analysis in 
Appendix D, getting a Gaussian approximation about each maximum. For large numbers 
of measurements these Gaussians act as delta-functions at the two angles ±$0 centered on 
<p = 0, as we have found in the experiment. Since there are two delta- functions then we find 
from Eq. (JTUJ) that 

Pm(+) = \[0- + c os(+$o)) + (1 + cos(-$ ))] = \0- + cos($ )) • (33) 
From Eq. (|72*j) in Appendix D, this is 

P m (+) = cos 2 ($ /2) = — ■ (34) 

771 

Of course, this is the same final value we got for P m {+) above in Eq. ()3())1 . We also comfirm 
the relation between a and $o : 



cos$ = sin a 
$o = 7r/2 - a 
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(35) 



as we found above by considering the CEP. As far as the general hidden direction (/?, 7) is 
concerned, all we can say is sin /3 cos 7 = cos$o with some ambiguity as to signs or 7r/2's in 
the angles. 

Let us see how this all works out in our experiment. In our first experiment we have 
cos 2 ($o/2) = 262/300 = 0.873 so that $ = 0.73. Thus the peaks of g will be at ±0.73 rad 
as we have found experimentally and is shown in Fig. El 

Can we can determine both the unknowns (3 and 7 from measurements and so remove 
the ambiguity between the two sides ±$ of the cone? The only possible way to do this is 
to use other measurement angles. If we continue to measure only in the transverse plane, 
however, it seems likely we will determine only a combination of (3 and 7 and not both; 
we would need to go out of the transverse plane to get both. We now consider changes in 
measurement angle <f)\. 



B. Changing the measurement angle: 

Suppose we now pick a 4> m +i — +&o, (we continue to take the original constant 0i to be 
zero), that is, we measure as close to the angle of one of the peaks as we can in the next 
measurement. We might miss it by 5. What this does is to eliminate completely (or almost) 
one of the peaks. We show this in Fig. 0] When working at a single angle we did not know 
on which side of the measurement direction the emerging polarization was. By moving in 
that direction (or away from it), we do determine the side. We can see this analytically by 
using the Gaussian approximation developed in Appendix D for the first m interations. If 
we move to +$0 then the value of g m +i($) is 

W ($) ~ [e-h^-^o? + e -W$+$o) 2 ] (1 + Vm+1 cos ($ _ $ + S )) 

« (1 + /7m+i cos 5) e -^-^ 2 + (1 + r, m+1 cos(2$ - 5)) e~^(*+*o) 2 (36) 

where 5 is our measurement inaccuracy. Now if it turned out that f] m +i = — 1 we will 
have most likely made the wrong choice (because if we are close to the correct angle, it is 
most likely that we will get an +1 result) and the wrong peak will be almost completely 
eliminated. If r) m+ i = +1, then the correct peak will be emphasized by 2 and the wrong 
peak given a smaller coefficient. As Fig.|3]shows, it works with just a few measurements (15 
in this case) at the new angle —0.73. We now know on which side of the azimuthal angle 
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lies. 



C. Iterations and Convergence: 



Next consider what continuing to measure m! more times at this new angle produces. 
If we have done a good guess at what the polarization angle is then we will get mostly up 
spin results and as factor containing two peaks from these last w! measurements (with very 
small separation 25) near $ . Then we will have 



e -|m(*-* ) 2 _|_ e -|m($+* ) 2 



e -i m '(<I>-$o+<5) 2 + e -|m'(*-$ -«) 2 



(37) 



We get a sharpening of the peak. (If 5 is too large we might still have two very closely-spaced 
peaks near $ .) 

Suppose we now have only the one peak after ra iterations and it is near $ = +$o- That 

is, 

g m {$) ~ e -!"»(*-*o) 2 (38) 

This result implies that the mean deviation in the measurements should go as sjl/m. We 
Fourier transform this Gaussian to find 



a™ = 2a r n e-^ cos $ ; = 2a™e~^ sin $ 



(39) 



so that 

tan$ m = tan$ ; A m = e~^. (40) 

and if we continue to make other measurements near this same angle $ = § + 5 m , where 
5 m is the error in our rath guess about the exact position of <&o, then 

1 



P(+) 



1 + e a™ (cos($ + 8 m ) cos $ + sin($ + 5 m ) sin $ ) 



1 + e 2 m (l + <5 m sin $ cos <3>o) 



1 + e 



1. 



(41) 



The equivalent simple-minded point of view is that the spin that was originally at polar angle 
a = $0 an d zero azimuthal angle is now is at polar angle a « 7r/2 and azimuthal angle 
+<3>o- But this is still the equivalent for purposes of measurement to a (half) cone of possible 
"real" 3D positions. 
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Eqs. ()38|) and (|40|) make explicit predictions of how the width of the peak and the ampli- 
tude will proceed with measurement number. In order to test those approximate forms let 
us set up the double peak in 10 measurements, produce a single peak by measuring near one 
peak angle for just 20 steps and then go back to the original measurement direction <p m = 
for a few hundred measurements. In starting anew we develop a different random phase, of 
course. The results are given in Figs. G3|7| The approximate equations work very well. 

For these last measurements, the Fourier transform of Eq. (}3T?j) is valid and using it in 
Eq. (JUJ) with <j) m = we find 



Pm( + ) = \ 



in 

2a™ 



I 1 + e a™ (cos $ ) . (42) 



Comparing with the more general point of view at unknown 3D emerging phase angle (/3, 7), 
given by Eq. (fTTj) we get 

sin/3 cos 7 = e~2^ cos $ (43) 

We might have (3 = $0 and 7 = 0or/3 = tt/2 and 7 = $0 but we can never "know" 
absolutely. (Indeed in conventional quantum mechanics the wave function is a mixture of all 
possible phase angles that give the same probability P(+).) Staying in the transverse plane 
never yields both the space angles, just this combination. We have now found the maximal 
information on the emerging phase angle that the present set of measurements can offer. 
Appendix A discusses how to extend the analysis to angles out of the transverse plane. 
However, we do not extend the simulation to that case here. 



V. CONCLUSIONS 



We have presented a simulated experiment where two Bose condensates, one with spin 
up and the other with spin down, have been allowed to interfere. The standard quantum 
mechanical view is that a Fock state \N + ,N_) is a linear combination of phase states. As 
we make a series of measurements of spin orientation we develop a state that is a mixture 
of such Fock states having varying numbers of up and down spins, but with constant total 
number of spins. As we proceed the state becomes a narrower mixture of phase states and 
the experiment gradually changes the state into nearly just one particular phase state. In 
convertional quantum mechanics, the measurement process is considered to have created the 
final phase value by continual wave function collapse. However, as one does the experiment, 
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one may have an equally appealing additional-variable view that the phase emerging was 
built-in all along and the experiment simply revealed what it was. For the experimentalist, 
one view works just as well as the other. The formalism developed in Ref. 2l|and generalized 
in Appendix A is particularly appropriate for showing this equivalence since it expresses the 
results of measurements on a Fock up-down state through an integral over the "conjugate" 
variable, the relative phase $ of the two states. The probability of any sequence of results 
then appears as a sum over the phase, exactly as in theories with additional variables. Each 
measurement produces a new initial state, with a new $ distribution; so that the change of 
this probability distribution can be obtained very easily for any sequence of measurements 
as we have done in our simulation. Thus we have been able to observe the emergence of 
the additional variable under the effect of successive quantum measurements. The Amherst 
experiment is found to be an ideal example of an actual experiment along these lines. 
LKB is UMR 8552 de TENS, du CNRS and de l'Universit Pierre et Marie Curie. 

VI. APPENDICES 

A. Extension to spin measurement at arbitrary directions 



In this section we generalize the calculations of Sec. II and Ref. |2JJ, where the directions 
of spin measurements were assumed to be in the transverse plane; we now assume that the 
directions are arbitrary. Instead of a single angle fa to characterize each measurement, we 
then need two angles and fa. (It is more convenient to come back to the standard notation 
and keep 9i for the polar angle; fa therefore now replaces the notation 6i of Ref. We 
also re-introduce orbital variables, which were ignored in Sec. II. We call \a) and \(3) the 
two spins states (up and down) and \& aij a(r) the corresponding field operators. The spin 
component along the direction of measurement is: 

a e M = cos9 [^(r)# a (r)-^(r)^(r)] + sinfl [e~^ tt+ (r)^(r) + e<* *J(r)tt a (r)" 

□ <44) 

and the projector [26J over the spin eigenstates: 



16 



with 97 = d=l and the local density is defined by: 

n(r) = *t (r)¥ (r) + ^(r)^(r). (46) 

The projector is then proportional to: 

P v {r;6, (p) ~ [1 + 77 cos 

+^sin^ [e - ** *t(r)*^(r)+e** *J(r)tf a (r)] . (47) 

As in §3.1 of Ref. 0, we assume that the system of spin particles is initially in the state: 

|^o) = \N a : u a ,a; N b : u b ,/3) (48) 

where the orbital states correspond to the wave functions (normalized to one): 

(r \u a ,b) = u aib (r). (49) 

The initial system is therefore simply the juxtaposition of a large number N a of particles 
condensed into the wave function u a (r) with a large number N b of particles condensed into 
the wave function u b (r). We now assume that M spin measurements are performed at points 
ri, r 2 , • • • ,r m in spin direction defined by angles (9\, 0i), (62, ^2), • • • , (0 m , <f> m ) and calculate 
the probability for obtaining a series of results T]i, r} 2 ,- ■ ■ ,rj m (all 77's are equal to ±1). The 
corresponding probabity is the average value in state l^o) of a product of P projectors 

The rest of the calculation is very similar to that Ref. |21( the only difference being the 
presence of the terms in rj cos 9 in Eq. ()47|) and the sin 9 factor, which do not change much the 
calculation. The same considerations apply on the conservation of the number of particles 
in each state and, in the limit where m <C N a , N b , we can express this conservation through 
an integral over a phase $. The probability then becomes proportional to the expression: 



/ — TTl I 1 +Vi co sOi] N a \u a (ri)\ 2 + [1 - rji cos 6i] N b \u b (ri)\ 2 

JO 2.7T TJy V 

+ rusine^NaN,, [^KNe^"^ + c.c.]|, (50) 

where c.c. stands for "complex conjugate." The second line ot this result can also be written 

as 



+ 277iSin0i y/N a N b \u a (ri) \ \u b (ri)\ cos [ffa) - fa - $] (51) 
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where £(r) is the relative phase of the two wave functions of the condensates: 



£(r) = arg [ii a (r)/u 6 (r)] . 



(52) 



This second line contains all the $ dependence of the probabilities; in other words it is the 
only origin of correlations between different measurements. It also contains the r dependence 
of the probability, which produces the fringes in space. 

The discussion of §§ 1.1 and 3.1 of Ref. shows that effect of the 2-th measurement 
on the $ distribution is merely to multiply this distribution by a function Di (<&) that is 
nothing but the content of the i-th curly brackect in Eq. (J5Uj) (with a normalization constant 
that is unimportant for our discussion here). This multiplication provides the evolution of 
the information on the relative phase $ that is obtained by this measurement. Whatever 
measurement parameters r$, (9i,fa) are arbitrarily chosen, the multiplying function Di (<&) 
is always the sum of a constant (first line of (|5()|l) plus a sinusoidal variation given by (|51jl. 
The former depends on r $ and 9i only, with the dependence involving only the densities of 
probabilities associated with the condensed states; the latter depends also on fa as well as 
the relative phase of the two states. Both depend in general on the result of the measurement 
r]i, as opposed to the situation for measurement in transverse directions only (9i = 7r/2). 

If, for instance, we assume that rji = +1, the maximum and the minimum of the contri- 
bution to the probability of the i-th measurement can be written as 



This provides the optimum value of 9i for each measurement position r if r\i = +1. But, 
if r]i = — 1, it is easy to see that cos#j/2 and sin 9i/2 are interchanged in (|53p. so that 
the the right hand side of (|H%j) now provides the inverse of tan#j/2; the optimum value of 
Qi corresponding to the two possible results are therefore symmetrical with respect to the 
horizontal plane. 

As a consequence, the choice of the optimum value of 0, is not easy in general, since it 
depends on the random result of an experiment that is not yet known when the apparatus 




(53) 



The best contrast will therefore be obtained if the minimum vanishes, that is if: 




(54) 
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is adjusted. If, nevertheless, one can locate the position of the measurement at a point 
where the two bosonic fields have the same intensity: 

N a \u a {v t )\ 2 = iVfcMr,)! 2 (55) 

the dilemma disappears: the optimum value of 0, is 7r/2, corresponding to a measurement 
performed in a direction of the transverse plane. In this case, the flexibility introduced by 
the new parameter 9i is useless. Nevertheless, condition Eq. (I55j) is not necessarily easy to 
meet, and may even be impossible for some configurations. Then, the optimization of the 
polar direction of measurement becomes relevant, but only if one already has a good idea in 
advance of what the most likely value of $ is, so that with an appropriate choice of (pi it is 
possible to infer what the result 77^ of the next measurement will be with a good probability. 
The conclusion is then that the flexibility introduced by Q{ may be useful, but only after 
a series a measurements has already been performed, so that <3> is already reasonably well 
known. 

B. Proof that <£ m = (mod ir) for measurement along a single direction 

Consider the situation in which the measurement angle remains the same in every 
measurement. Then we can Fourier transform 

m— 1 

= H(l + ViCOs(® -0)). (56) 
1=1 

in a series in cos(0 — $) rather than in cos<3? and sin$. That is 

g m (§) = a ™ + J2c™ cos(g($ - 0)) = a™ + ^ (c™ cos $ cos + c™ sin $ sin 0) . (57) 
from which we find 

/>„(+)- 1 + ^. (58) 
By comparison with Eq. (|13j) we see that 

a™ = c ™ cos 0; b™ = c™ sin 0, (59) 

which, by Eq. (jl5j) . shows that tan$ m = tan0 or 

$ m = or + 7T, (60) 
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and that c™ /2a™ = sin a m cos($ m — 0) = ± sin a m . This also gives 

P m {+) = — = ~(l±sina m ). (61) 
m 2 

where the + sign occurs with m + > m_ and the — sign in the opposite case so that sin a m 
is positive. 



C. Analytic results for constant measurement angle 

When all measurement angles are the same value we can with full generality take that 
angle equal to zero and write 

m—l 

g m m = n( i +^ cos ( $ )) = ( i + cos ( $ )) m+ (i-cos($)r- 

i=l 

/ <f>\ m + / $\ m - 

= 2 m (^cos-j (sin -J (62) 
Then the Fourier transforms are found from integral tables j3] to be 

<C = — / d® cos — sin — 



2tx h V 2 J \ 2 



2 - r (m+ + 1) r ( 



/»- + ! 



7r r (m + 1) 



a™ = — y d$(cos— J ( v sin "^"' Cos * 



(63) 



2 m+i r (m+ + 1) r (m_ + 1) - r (m+ + i) r (m_ + 1 



(64) 



7i r (m + 2) 

with 6™ = by symmetry. Taking the ration of these and using the definition of sin a m gives 

W?\ \m+-m-\ , . 

2a™ m + 1 

This result allows us to understand the behavior of A m as seen in Fig. 1. We notice that 
the plot has small upward curves followed by abrupt downward jumps. We find that the 
upward sweep is a sequence of all rj = 1 results while a downward jump is a single rj = — 1. 
Suppose, as occurs in the figure that m + is considerably larger than m„. Then when a new 
up spin result occurs the change in A m is given by 

. m + + 1 — m_ m + — m_ 2m_ + 1 

A + = = (66) 

m + 2 m+1 (to + 2) (to + 1) v ' 
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while, if a down spin occurs in the measurement, the change in this quantity is 

= m + - m_ - 1 _ m + - m_ = 2m + + 1 

^~ m + 2 m + 1 ~~ (m + 2)(m+l) 

The second quantity, the jump down, is much greater in magnitude than the previous upward 
move, giving the peculiar shape of the curve. 



D. The peaks in g(<&) for the spin phase 

Here we analyze in the case of all equal measurement axes fa = to show it has 
two sharp peaks. We can take the single measurement angle to be zero with generality (or 
work with $' = $ — fa). The location of these peaks corresponds to what we have found 
numerically. Write g(&) in the half- angle form 

g{S) = (cos 2 [($ )/2])" i+ (sin 2 [($ )/2]) m - . (68) 

Set the derivative of the logarithm of this to zero to determine the position of the maxima. 

d$ cos[($)/2] sin[($)/2] v ; 

which is equivalent to 

cos 2 [($)/2l = (70) 
m 

This equation has two solutions in the range < $ < 2n. These are at 

$ = ±$o (71) 

where <3> is in the range [0, n] or 

cos 2 [$ /2] = (72) 
m 

Thus the two peaks occur symmetrically about the value of $ = just as we have found 
"experimentally." Correspondingly we have 

sin 2 [$ /2] = — . (73) 
m 

We get a Gaussian approximation to g by Taylor expanding lng about ±$o- We use 

d 2 lng{$) 



d<S> 2 
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= -m, (74) 

*0 



so that 



□ 



9(*) 



m + \ m + /m_ 
m J \ m 



e -im($-$ ) 2 _|_ e -|m($+$ ) 2 



(75) 



For large m these act like delta-functions. 
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FIG. 1: The cone of equal probability. Any spin at angles /3, 7 on this cone will give the same 
probability of being up along the measurement axis, taken as x in this figure. 
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FIG. 2: The amplitude A = sin a as a function of interation step when all measurements have been 
made at a single angle. The final asymptotic value is random. 
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FIG. 3: The angular distribution g($) as a function of angle for three iteration step lengths, 10 
steps (dashed line), 150 steps (dottted line), and 300 steps (solid line). For a single measuring 
angle this always has two equal peaks corresponding to the intersection of the spin cone with the 
transverse plane. The peaks narrow with step length. 
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FIG. 4: Typical g(&) after the second peak has been eliminated by a few measurements at the 
positive peak. 
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FIG. 5: The amplitude A m as a function of interation step when all measurements, after eliminating 
the second peak of g($>) early on, are back at the original measurement angle (f>\ = 0. The dotted 
line is the approximation e _1 / 2m . 
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FIG. 6: The azimuthal angle <£ m « $o as a function of interation step when all measurements, 
after eliminating the second peak of g(&) early on, are back at the original measurement angle 
01 = 0. 
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FIG. 7: The mean deviation (width of g(3>)) as a function of interation step when all measurements, 
after eliminating the second peak of g(&) early on, are back at the original measurement angle 
<pi = 0. The dotted line is the approximate form 
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